Search for: All records

ABSTRACT JWST has revealed a large population of UV-bright galaxies at $$z\gtrsim 10$$ and possibly overly massive galaxies at $$z\gtrsim 7$$, challenging standard galaxy formation models in the ΛCDM cosmology. We use an empirical galaxy formation model to explore the potential of alleviating these tensions through an Early Dark Energy (EDE) model, originally proposed to solve the Hubble tension. Our benchmark model demonstrates excellent agreement with the UV luminosity functions (UVLFs) at $$4\lesssim z \lesssim 10$$ in both ΛCDM and EDE cosmologies. In the EDE cosmology, the UVLF measurements at $$z\simeq 12$$ based on spectroscopically confirmed galaxies (eight galaxies at $$z\simeq 11\!-\!13.5$$) exhibit no tension with the benchmark model. Photometric constraints at $$12 \lesssim z\lesssim 16$$ can be fully explained within EDE via either moderately increased star-formation efficiencies ($$\epsilon _{\ast}\sim 3\!-\!10\ \hbox{per cent}$$ at $$M_{\rm halo}\sim 10^{10.5}{\, \rm M_\odot }$$) or enhanced UV variabilities ($$\sigma _{\rm UV}\sim 0.8\!-\!1.3$$ mag at $$M_{\rm halo}\sim 10^{10.5}{\, \rm M_\odot }$$) that are within the scatter of hydrodynamical simulation predictions. A similar agreement is difficult to achieve in $$\Lambda$$CDM, especially at $$z\gtrsim 14$$, where the required $$\sigma _{\rm UV}$$ exceeds the maximum value seen in simulations. Furthermore, the implausibly large cosmic stellar mass densities inferred from some JWST observations are no longer in tension with cosmology when the EDE is considered. Our findings highlight EDE as an intriguing unified solution to a fundamental problem in cosmology and the recent tensions raised by JWST observations. Data at the highest redshifts reached by JWST will be crucial for differentiating modified galaxy formation physics from new cosmological physics. 

Abstract We analyze the first cosmological baryonic zoom-in simulations of galaxies in dissipative self-interacting dark matter (dSIDM). The simulations utilize the FIRE-2 galaxy formation physics with the inclusion of dissipative dark matter self-interactions modeled as a constant fractional energy dissipation (f_diss= 0.75). In this paper, we examine the properties of dwarf galaxies withM_*∼ 10⁵–10⁹M_⊙in both isolation and within Milky Way–mass hosts. For isolated dwarfs, we find more compact galaxy sizes and promotion of disk formation in dSIDM with (σ/m) ≤ 1 cm²g⁻¹. On the contrary, models with (σ/m) = 10 cm²g⁻¹produce puffier stellar distributions that are in tension with the observed size–mass relation. In addition, owing to the steeper central density profiles, the subkiloparsec circular velocities of isolated dwarfs when (σ/m) ≥ 0.1 cm²g⁻¹are enhanced by about a factor of 2, which are still consistent with the kinematic measurements of Local Group dwarfs but in tension with the Hirotation curves of more massive field dwarfs. Meanwhile, for satellites of Milky Way–mass hosts, the median circular velocity profiles are marginally affected by dSIDM physics, but dSIDM may help promote the structural diversity of dwarf satellites. The number of satellites is slightly enhanced in dSIDM, but the differences are small compared with the large host-to-host variations. In conclusion, the dSIDM models with (σ/m) ≳ 0.1 cm²g⁻¹,f_diss= 0.75 are in tension in massive dwarfs (M_halo∼ 10¹¹M_⊙) due to circular velocity constraints. However, models with lower effective cross sections (at this halo mass/velocity scale) are still viable and can produce nontrivial observable signatures. 

ABSTRACT Recent observations with JWST have uncovered unexpectedly high cosmic star formation activity in the early Universe, mere hundreds of millions of years after the big bang. These observations are often understood to reflect an evolutionary shift in star formation efficiency (SFE) caused by changing galactic conditions during these early epochs. We present FIREbox$$^{\it HR}$$, a high-resolution, cosmological hydrodynamical simulation from the Feedback in Realistic Environments (FIRE) project, which offers insights into the SFE of galaxies during the first billion years of cosmic time. FIREbox$$^{\it HR}$$ re-simulates the cosmic volume ($L=22.1$ cMpc) of the original FIREbox run with eight times higher mass resolution ($$m_{\rm b}\sim {}7800\, M_\odot$$), but with identical physics, down to $$z\sim {}6$$. FIREbox$$^{\it HR}$$ predicts ultraviolet (UV) luminosity functions in good agreement with available observational data. The simulation also successfully reproduces the observed cosmic UV luminosity density at $$z\sim {}6{\!-\!}14$$, demonstrating that relatively high star formation activity in the early Universe is a natural outcome of the baryonic processes encoded in the FIRE-2 model. According to FIREbox$$^{\it HR}$$, the SFE–halo mass relation for intermediate mass haloes ($$M_{\rm halo}\sim {}10^9{\!-\!}10^{11}\, {\rm M}_\odot$$) does not significantly evolve with redshift and is only weakly mass-dependent. These properties of the SFE–halo mass relation lead to a larger contribution from lower mass haloes at higher z, driving the gradual evolution of the observed cosmic UV luminosity density. A theoretical model based on the SFE–halo mass relation inferred from FIREbox$$^{\it HR}$$ allows us to explore implications for galaxy evolution. Future observations of UV faint galaxies at $$z\gt 12$$ will provide an opportunity to further test these predictions and deepen our understanding of star formation during Cosmic Dawn. 

Abstract Using cosmological hydrodynamical zoom-in simulations, we explore the properties of subhalos in Milky Way analogs that contain a subcomponent of atomic dark matter (ADM). ADM differs from cold dark matter (CDM) due to the presence of self-interactions that lead to energy dissipation, analogous to standard model baryons. This model can arise in dark sectors that are natural and theoretically motivated extensions to the standard model. The simulations used in this work were carried out usingGIZMOand utilize the FIRE-2 galaxy formation physics in the standard model baryonic sector. For the parameter points we consider, the ADM gas cools efficiently, allowing it to collapse to the center of subhalos. This increases a subhalo’s central density and affects its orbit, with more subhalos surviving small pericentric passages. The subset of subhalos that host satellite galaxies have cuspier density profiles and smaller stellar half-mass radii relative to CDM. The entire population of dwarf galaxies produced in the ADM simulations is more compact than those seen in CDM simulations, unable to reproduce the entire diversity of observed dwarf galaxy structures. Additionally, we also identify a population of highly compact subhalos that consist nearly entirely of ADM and form in the central region of the host, where they can leave distinctive imprints in the baryonic disk. This work presents the first detailed exploration of subhalo properties in a strongly dissipative dark matter scenario, providing intuition for how other regions of ADM parameter space, as well as other dark sector models, would impact galactic-scale observables. 

The distribution of offsets between the brightest cluster galaxies of galaxy clusters and the centroid of their dark matter distributions is a promising probe of the underlying dark matter physics. In particular, since this distribution is sensitive to the shape of the potential in galaxy cluster cores, it constitutes a test of dark matter self-interaction on the largest mass scales in the universe. We examine these offsets in three suites of modern cosmological simulations; IllustrisTNG, MillenniumTNG and BAHAMAS. For clusters above , we examine the dependence of the offset distribution on gravitational softening length, the method used to identify centroids, redshift, mass, baryonic physics, and establish the stability of our results with respect to various nuisance parameter choices. We find that offsets are overwhelmingly measured to be smaller than the minimum converged length scale in each simulation, with a median offset of in the highest resolution simulation considered, TNG300-1, which uses a gravitational softening length of . We also find that centroids identified via source extraction on smoothed dark matter and stellar particle data are consistent with the potential minimum, but that observationally relevant methods sensitive to cluster strong gravitational lensing scales, or those using the the “light traces mass” approach, in this context meaning gas is used as a tracer for the potential, can overestimate offsets by factors of $\sim 10$and $\sim 30$, respectively. This has the potential to reduce tensions with existing offset measurements which have served as evidence for a nonzero dark matter self-interaction cross section. 

ABSTRACT Both observations and simulations have shown strong evidence for highly time-variable star formation in low-mass and/or high-redshift galaxies, which has important observational implications because high-redshift galaxy samples are rest-ultraviolet (rest-UV) selected and therefore particularly sensitive to the recent star formation. Using a suite of cosmological ‘zoom-in’ simulations at z > 5 from the Feedback in Realistic Environments project, we examine the implications of bursty star formation histories for observations of high-redshift galaxies with JWST. We characterize how the galaxy observability depends on the star formation history. We also investigate selection effects due to bursty star formation on the physical properties measured, such as the gas fraction, specific star formation rate, and metallicity. We find the observability to be highly time-dependent for galaxies near the survey’s limiting flux due to the star formation rate variability: as the star formation rate fluctuates, the same galaxy oscillates in and out of the observable sample. The observable fraction $$f_\mathrm{obs} = 50~{{\ \rm per\ cent}}$$ at z ∼ 7 and M⋆ ∼ 108.5–$$10^{9}\, {\rm M}_{\odot }$$ for a JWST/NIRCam survey reaching a limiting magnitude of $$m^\mathrm{lim}_\mathrm{AB} \sim 29{\!-\!}30$$, representative of surveys such as JADES and CEERS. JWST-detectable galaxies near the survey limit tend to have properties characteristic of galaxies in the bursty phase: on average, they show approximately 2.5 times higher cold, dense gas fractions and 20 times higher specific star formation rates at a given stellar mass than galaxies below the rest-UV detection threshold. Our study represents a first step in quantifying selection effects and the associated biases due to bursty star formation in studying high-redshift galaxy properties. 

ABSTRACT JWST observations have revealed a population of galaxies bright enough that potentially challenge standard galaxy formation models in the Λ cold dark matter (ΛCDM) cosmology. Using a minimal empirical framework, we investigate the influence of variability on the rest-frame ultra-violet (UV) luminosity function of galaxies at z ≥ 9. Our study differentiates between the median UV radiation yield and the variability of UV luminosities of galaxies at a fixed dark matter halo mass. We primarily focus on the latter effect, which depends on halo assembly and galaxy formation processes and can significantly increase the abundance of UV-bright galaxies due to the upscatter of galaxies in lower-mass haloes. We find that a relatively low level of variability, σUV ≈ 0.75 mag, matches the observational constraints at z ≈ 9. However, increasingly larger σUV is necessary when moving to higher redshifts, reaching $$\sigma _{\rm UV} \approx 2.0\, (2.5)\, {\rm mag}$$ at z ≈ 12 (16). This implied variability is consistent with expectations of physical processes in high-redshift galaxies such as bursty star formation and dust clearance during strong feedback cycles. Photometric constraints from JWST at z ≳ 9 therefore can be reconciled with a standard ΛCDM-based galaxy formation model calibrated at lower redshifts without the need for adjustments to the median UV radiation yield. 

Abstract Recent discoveries of a significant population of bright galaxies at cosmic dawn $\left(z\gtrsim 10\right)$have enabled critical tests of cosmological galaxy formation models. In particular, the bright end of the galaxys’ UV luminosity functions (UVLFs) appear higher than predicted by many models. Using approximately 25,000 galaxy snapshots at 8 ≤z≤ 12 in a suite of FIRE-2 cosmological “zoom-in” simulations from the Feedback in Realistic Environments (FIRE) project, we show that the observed abundance of UV-bright galaxies at cosmic dawn is reproduced in these simulations with a multichannel implementation of standard stellar feedback processes, without any fine-tuning. Notably, we find no need to invoke previously suggested modifications, such as a nonstandard cosmology, a top-heavy stellar initial mass function, or a strongly enhanced star formation efficiency. We contrast the UVLFs predicted by bursty star formation in these original simulations to those derived from star formation histories (SFHs) smoothed over prescribed timescales (e.g., 100 Myr). The comparison demonstrates that the strongly time-variable SFHs predicted by the FIRE simulations play a key role in correctly reproducing the observed, bright-end UVLFs at cosmic dawn: the bursty SFHs induce order-or-magnitude changes in the abundance of UV-bright (M_UV≲ −20) galaxies atz≳ 10. The predicted bright-end UVLFs are consistent with both the spectroscopically confirmed population and the photometrically selected candidates. We also find good agreement between the predicted and observationally inferred integrated UV luminosity densities, which evolve more weakly with redshift in FIRE than suggested by some other models. 

ABSTRACT Using high-resolution cosmological radiation-hydrodynamic (RHD) simulations (thesan-hr), we explore the impact of alternative dark matter (altDM) models on galaxies during the Epoch of Reionization. The simulations adopt the IllustrisTNG galaxy formation model. We focus on altDM models that exhibit small-scale suppression of the matter power spectrum, namely warm dark matter (WDM), fuzzy dark matter (FDM), and interacting dark matter (IDM) with strong dark acoustic oscillations (sDAO). In altDM scenarios, both the halo mass functions and the ultraviolet luminosity functions at z ≳ 6 are suppressed at the low-mass/faint end, leading to delayed global star formation and reionization histories. However, strong non-linear effects enable altDM models to ‘catch up’ with cold dark matter (CDM) in terms of star formation and reionization. The specific star formation rates are enhanced in halos below the half-power mass in altDM models. This enhancement coincides with increased gas abundance, reduced gas depletion times, more compact galaxy sizes, and steeper metallicity gradients at the outskirts of the galaxies. These changes in galaxy properties can help disentangle altDM signatures from a range of astrophysical uncertainties. Meanwhile, it is the first time that altDM models have been studied in RHD simulations of galaxy formation. We uncover significant systematic uncertainties in reionization assumptions on the faint-end luminosity function. This underscores the necessity of accurately modeling the small-scale morphology of reionization in making predictions for the low-mass galaxy population. Upcoming James Webb Space Telescope imaging surveys of deep lensed fields hold potential for uncovering the faint low-mass galaxy population, which could provide constraints on altDM models. 

ABSTRACT Cosmological simulations serve as invaluable tools for understanding the Universe. However, the technical complexity and substantial computational resources required to generate such simulations often limit their accessibility within the broader research community. Notable exceptions exist, but most are not suited for simultaneously studying the physics of galaxy formation and cosmic reionization during the first billion years of cosmic history. This is especially relevant now that a fleet of advanced observatories (e.g. James Webb Space Telescope, Nancy Grace Roman Space Telescope, SPHEREx, ELT, SKA) will soon provide an holistic picture of this defining epoch. To bridge this gap, we publicly release all simulation outputs and post-processing products generated within the thesan simulation project at www.thesan-project.com. This project focuses on the z ≥ 5.5 Universe, combining a radiation-hydrodynamics solver (arepo-rt), a well-tested galaxy formation model (IllustrisTNG) and cosmic dust physics to provide a comprehensive view of the Epoch of Reionization. The thesan suite includes 16 distinct simulations, each varying in volume, resolution, and underlying physical models. This paper outlines the unique features of these new simulations, the production and detailed format of the wide range of derived data products, and the process for data retrieval. Finally, as a case study, we compare our simulation data with a number of recent observations from the James Webb Space Telescope, affirming the accuracy and applicability of thesan. The examples also serve as prototypes for how to utilize the released data set to perform comparisons between predictions and observations.